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Abstract 

X-ray binary Cygnus X-l is microquasar containing historically first candidate to black hole. 
Paradigm of magnetic disc accretion dominates in theoretical models describing processes taking 
place in objects containing black holes such as microquasars and active galactic nuclei. Nevertheless 
up to now there were no reliable measurements of magnetic fields in these systems. The first prediction 
of the Cyg X-l magnetic field was done by V.F.Shvartsman . He wrote that Cyg X-l X-ray emission 
millisecond flickering evidences the presence of a black hole and points to the magnetic field role in 
accretion onto a black hole 2,3 . From that times there were many of attempts to search for the Cyg 
X-l magnetic field but all of these efforts indicated upper limits only. Our VLT FORS1 2007 and 
2008 observations revealed a presence of a magnetic field in the system. For the first time we obtained 
on the level of 6 standard deviations (a) a magnetic field of the order of 130 G on the surface of the 
Cyg X-l optical component (O-supergiant) and observational estimation on a 4<r level of a magnetic 
field on a outer part of accretion structure (about 600 G) in accordance with theoretical prediction. 
Scaling this field value to black hole vicinity we showed that the field is strong enough to explain the 
X-ray millisecond flickering and can explain it. Our result presents the first direct determination of 
magnetic field in accreting disk around a black hole. 

Cyg X-1/HDE226868 is a X-ray binary system with the orbital period P=5.6 d , whose relativistic com- 
ponent is the first candidate black hole (BH). The optical component (09.7 lab supergiant) is responsible 
for about 95% of the system optical luminosity. The remaining 5% are due to the accretion structure 
(disc and surrounding gas) near BH. 

Though the investigation of Cyg X-l are being carried out over 40 years and have resulted in ~1000 
publications, geometrical and physical parameters the system remain unclear. The same is also true for 
some phenomena observed in this system including e.g. long time periodic and aperiodic variations and 
flares. 

Both linear and circular polarizations was detected in Cyg X-l optical continuum and investigated 
in 1970s. Linear interstellar and circumstellar polarization reaches a value of ~5% 4 . Its strength and 
position angle change on the scale of years 5 . Kemp 6 detected a component of an amplitude ~0.25% 
which is variable over the orbital phase 4>. This component shows complicated and variable dependence 
from 4>. Similar behaviour was found for circular polarization discovered in 1972 7,8 . Interstellar circular 
polarization does not exceed 0.04%, and variable component with the period of 2.8/5.6 d is about 0.02% 
9 . While intrinsic circular polarization is most probably generated by a magnetic field, the intrinsic linear 
polarization in Cyg X-l is usually explained by electron scattering on non-symmetrical gas structures 
(e.g. 10 " 13 ). 

Theoretical estimations of the strength of the magnetic field in Cyg X-l were based on optical polar- 
ization 14 ' 15 . Upper limits B < 350 G for the optical component and B < 500 G for the outer part of 
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Table 1: Magnetic fields from VLT spectropolarimetrical observations of X-ray binary Cyg X-l. JD - 
Julian date of the middles of the object observations; Orbital phases </>, according cphcmcrcdcs 22 , = 
- optical component is ahead; a - standard deviation, \(B z )/a\ - significance. 
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accretion disc were found. The 6-m telescope spectropolarimetric observations show B < 1000 G for the 
outer part of the disc 16 . 

Our method of the study of the presence of a magnetic field is based on measurements of the circular 
polarization (the Stokes parameter V) in optical spectra produced by Zeeman effect. The method of the 
determination of the mean (averaged over the picture plane) longitudinal magnetic field (B z ) is described 
in full detail by 17 . The method is statistical: to increase the sensibility there are used simultaneously all 
observed spectral lines. Value (B z ) is obtained 17-19 from the slope of least squares linear regression of 

Y = V/I 

and 

X = -4.67 x 10~ 13 x g eff x (B z ) x A 2 x (1//) x dl/dX. 

Here I is intensity of radiation, A is wavelength in angstroms, t/ e // is effective Lande factor. 

Our spectropolarimetrical observations were conducted with the European South Observatory Very 
Large Telescope (VLT) 8.2 m (Cerro Paranal, Chile) in service mode with the FORS1 spectrograph in 
the range 3680-5129 A, spectral resolution R=4000, signal-to-noise ratio S/N =1500 - 3500 (for spectra 
of intensity) in 2007 from June 18 to July 9 and in July, 2008 (see Table 1). The system Cyg X-l was at 
that time in its X-ray "hard state". 13 spectropolarimetric spectra with exposure time of ~1 hour were 
obtained during 13 nights. For our observations we adopt effective Lande-factor g e //=1.07 according 20 . 

The used method has been already applied (and carefully tested) in previous studies of bright magnetic 
stars 17-19 which usually do not have significant interstellar or intrinsic linear polarization and have rather 
strong (B z ) =500-2000 G. In contrast, Cyg X-l has (B z ) weaker and strong interstellar / circumstellar 
linear polarization. For this reason we had to meet some precautions in the magnetic field measurements 
of this system and to adapt the method for such conditions. 

Before (B z ) calculations we clean carefully V/I spectra from any features which could distort results. 
We excluded features alien to photosphere of Cyg X-l optical component: 1) the wavelengths of inter- 
stellar lines; 2) defects (including weak residual cosmic ray tracks remained after standard observation 
processing); 3) Hell 4686 emission line and 4) emission components of the lines with strong P Cyg effect. 
We did not found the pollution by telluric lines in our spectra. 

The continua of the recorded V/I spectra show slight slopes. The continuum level drop 0.05% - 0.15% 
within our spectral range 3680-5129 A varied from night to night. This behaviour cannot be explained by 
Cyg X-l interstellar or/and intrinsic circular polarization. The probable reason could be a presence of a 
cross-talk between linear and circular polarization within the FORS1 analysing equipment. In agreement 



2 



with previous studies of Cyg X-l linear polarization in optical range 5 this effect is observed only in the 
continuum and does not distort S-shape V profiles of spectral lines caused by Zeeman effect (hereinafter 
Zeeman S-waves). 

To avoid any impact of the continuum slope on our (B z ) measurements, we subtracted linear trends 
from V/I spectra. After removal of slopes the (B z ) becomes lower by 20-80 G depending on the slope 
values (all (B z ) corrections are negative). 

We normalized I-spectra by pseudo-continuum following 18 > 19 . I-continuum is produced by the source 
energy distribution, interstellar reddening, broad diffuse interstellar bands (DIBs) as well as atmospheric 
extinction and detector sensitivity. Its slope reaches \d(log(I(X))/d(log(X))\ ~ 20. The slope removing 
gives (B z ) correction up to ~ 20 G. It is usually less than the statistical errors a((B z )) ^ 20 — 30 G, see 
table 1. 

After abovementioned reductions the residual deviations of least squares linear regression become to 
follow to the Gauss function satisfactorily up to ±3.6 tr(V/I), where <t(V/I) is standard deviation of V/I. 
It is mean that the level of significance corresponds to Gauss statistics now. Small number of points has 
larger deviations. Weak cosmic ray tracks can be a source of such deviations. Therefore we excluded 
from our analysis any pixel showing residual deviations exceeding 3.6 x a (V/I). 

The results of our measurements are presented in table 1 and figure 1 . To verify our results we used 
some tests: 1) each spectrum was divided in two halves at mid-wavelength; we checked that (Z? z )-values 
determined over each half separately were in agreement within error-bars. 2) We repeated (B z ) calculation 
using fragments of spectra, which include strong absorption lines (deeper than 4%) only; ~l/3 spectral 
points were used; (B z )- values were consistent with the previous measurements using the whole spectral 
regions within 1.5 sigma ((B z )). 3) Zeeman S-waves for the strongest lines were found. In Fig. 2a we 
show an example of a distinct Zeeman feature in the Hel line at A4026A. 

We should note that the element overabundance on the factor from 2 to 10 in the Cyg X-l optical 
component stellar atmosphere 22_23 enforces the spectral lines and Zeeman S-waves. It increases accuracy 
of (B z ) measurements. 

As a next step of our study we investigated the spectral line HeIIA4686A separately. Due to the 
presence of a strong emission component in the line profile it was omitted from the earlier analysis. In 
fact, this line has compound profile consisting of absorption (originating in the stellar photosphere) and 
emission (originating in the accretion structure) components. Certainly, the accuracy of the measurement 
of the magnetic field using just one line is considerably worse compared with measurement using the 
whole spectrum. Nevertheless, our analysis shows that for two spectra accuracies of estimations exceed 
4cr level: (5^)=-780+/-177 G in 2007 for the orbital phase <^=0.65 and (B 2 )=449+/-112 G in 2008 for 
</?=0.43. Zeeman S-wave in V-spectrum smoothed over 3 A and its correspondence to the dI(A)/dA wave 
is presented on figure 2b. 

To find He II 4686 A line formation regions we constructed Doppler tomogram (the binary system 
image in velocity space) on the base of our VLT observational data. We used new Doppler tomogram 
reconstruction technique worked out by Agafonov , so-called Radioastronomical Approach (RA). This 
RA method includes an effective CLEAN procedure and allows to reconstruct well the 2D velocity field 
with very small number of ID profiles (5-10 spectra may be sufficient), see as example the figure 1 in 25 . 
The tomography map constructed on the base of all 13 HeIIA4686 A line VLT profiles obtained by us in 
2007 and 2008 is presented in figure 3. It shows that HcIIA4686 A line emission regions are located near 
the point LI in the Roche lobe model and near the "hot spot" or "hot line" 26 on the outer part of the 
accretion structure. But for the different observational seasons Cyg X-l tomogram maps may differ from 
one another - matter flow changes at the scale of years 27-28 . Therefore we constructed also tomograms 
using our VLT-profilcs separately for 2007 and 2008. They show a similar result. Consequently (B z ) 
derived from Hell emission line is located in outer parts of accretion structure. Its values 400-800 G are 
in agreement with estimation 15 . 

Our main conclusions are the following: We discovered a longitudinal magnetic field averaged over the 
picture plane (B z ) ~ 100 G in the photosphere of Cyg^X-1 optical component. Real magnetic field can 
exceed (B z ). Magnetic field was detected at high confidence level (about 6a) near orbital phase (p—0A 
in 2007, and </?=0.5 in 2008. Dependence of (B z ) from ip is more complicated than for magnetic dipolc 
and is probably changed during one year (figure 1). The magnetic field structure variation may be the 
reason of some long-time variations of matter flowing process in this binary system. 

Further we found (B z ) ~ 600 G in outer parts of the accretion structure surrounding the BH. It is in 
agreement with Shvartsman's ideas 3 , that gas stream carries the magnetic field to the accretion structure 
and the gas is compressed by a factor of ~1 due to the interaction with the structure of the outer rim. 
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Gas density is increased and magnetic field is increased up to B ~ 600 G at a distance from BH 6 x 10 
cm = 2 x 10 5 x R g 29 (R g is gravitation radius). According to Shakura-Sunyaev 30 magnetized accretion 
disc standard model at 3 x R g B ~ 10 9 G. Taking into account radiative pressure predominance inside 
~10 — 20 Kg, we get B(3 x R g ) ~(2 — 3)xl0 8 G. The measured value of the magnetic field strength at the 
marginal orbit of the Cyg X-l black hole corresponds quite well to the Magnetic Coupling model with 
equipartition between kinetic and magnetic energy densities. 

If the X-ray millisecond flickering is related to the magnetic nature, then the accreting matter magnetic 
energy flux must exceed the X-ray emission fluctuating component luminosity. X-ray emission originates 
at R < 30 x R g . Inside the sphere of this radius the magnetic energy amounts to 10 40 erg. The radial 
velocity of magnetized plasma at 30 x R g in Shakura-Sunyaev accretion disk is ~1.5 km/s (we adopt 
viscosity parameter a — 1, because magnetic viscosity is big 31 ~ 32 ). The time of matter fall is ^1000 s. 
Magnetic energy flux is 10 37 erg/s which is equal or exceed the flickering component power (0.5 — 1) x 10 37 
erg/s. So, magnetic energy dissipation permits to account for the X-ray flickering. 

From abovementioned estimations the Cyg X-l black hole magnetic moment is about 10 30 G x cm 3 . 
According to 33 such object belongs to Magnetic Extremely Compact Object (MECO) class. 

So VLT FORS1 2007-2008 observations permit to detect the presence of a magnetic field in Cyg X-l. 
It is the pioneer measurement in black hole systems. The field can be responsible for X-ray millisecond 
flickering. Our result points to necessity of taking into account of magnetic field impact on the matter 
fluid structure in Cyg X-l. 
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Figure captions 

Figure 1. The mean longitudinal magnetic field of the Cyg X-l optical component (B z ) (in 
Gauss) vs. the orbital phase <f>. Full squares: 2007 data; full triangles: 2008 data. Error bars are 
68% confidence intervals. The phases of orbital period 5.6 days were calculated with ephemeredes 21 , <f> 
— corresponds to the time, when BH component is located behind the optical component (09.7 lab 
supergiant). a, The graph for 2007. b, The graph for 2008. c, The graph for 2007 and 2008 together. 

Figure 2. Examples of Zeeman S-waves. Two examples showing accordance of wavelength depen- 
dence of intensity normalized to continuum I/I c to Zeeman S-waves of V/I spectra are given, a, An 
example for absorption spectral lines originated in atmosphere of O-star component of Cyg X-l binary 
system: Zeeman S-wave (down box) in the region of Hel A4026A line (upper box) observed July 9, 2007 
when (B z ) was high (see table 1). b, Evidence of magnetic field in accreting gas following from Hell 
A4686A emission line analyzes for June 18, 2007. The line on the upper panel shows the observed Hell 
A4686A spectral line profile I/I c . Horizontal line is pseudo-continuum level. Solid line on the down 
panel shows V/I spectrum smoothed over 3 A. The dashed line on the panel shows in arbitrary units the 
expected Zeeman S-wave shape (dI/dA)/I (V/I ~ (dI/dA)/I), where I was smoothed over 3 A. Vertical 
solid line shows the centre of Hell A4686 A emission. 

Figure 3. Doppler tomogram of Cyg X-l in Hell A4686 A. It is constructed on the base of 13 
VLT spectra from 2007-2008 and shows brightness of emission in Hell A4686 A spectral line on velocities 
plane (V x ,V y ). Thin lines arc isophots of different levels shown on the right part of figure (in relative 
units). The zero level corresponds to "pseudo-continuum"; negative values correspond to absorption and 
positive ones - to emission. The dashed line shows the Roche lobe of black hole (BH) . The almost filling 
its Roche lobe optical component is drown by solid line. The Roche lobes are constructed for the mass 
ratio q = Mx/Mo = 1/3. Here Mx is the mass of X-ray emitting component (BH) and Mo is the mass 
of optical component (O-star). Ovals represent outer parts of the accretion disks with radii = 0.2 and 
0.25 of the distance between mass centers of the components. 
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